Investigation of a Toy Model for Frustration in Abelian Lattice Gauge Theory 
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We introduce a lattice model with local U(l) gauge sym- 
metry which incorporates explicit frustration in d > 2. The 
form of the action is inspired from the loop expansion of the 
fermionic determinant in standard lattice QED. We study 
through numerical simulations the phase diagram of the 
model, revealing the existence of a frustrated (antiferromag- 
netic) phase for d — 3 and d — 4, once an appropriate order 
parameter is identified. 



I. INTRODUCTION 

The importance of frustration and disorder is well 
known to condensed matter physicists working in the 
field of spin-glasses and related systems [Q . In these sys- 
tems one can find a variety of unusual phases which, in 
some cases, present completely new phenomena as, for 
instance, the existence of a multiplicity of vacua not re- 
lated by a symmetry of the action. It is interesting to 
ask whether similar phenomena could happen in a field 
theory and what will be the physical consequences. 

We address here the issue of frustration in lattice gauge 
theories. The motivation for this work relies on the fact 
that, as well known, realistic gauge theories with dynam- 
ical fermions, in contrast to gauge theories with funda- 
mental scalar fields, are characterized by a frustrated 
effective action. Frustration here must be understood 
in the most naive way, i.e., after the integration of the 
Grassmann degrees of freedom, the effective fermion ac- 
tion has contributions which can be seen as competing 
interactions. 

Even if spin-glass like phases in gauge theories have 
been found time ago we are interested, as stated before, 
in more realistic gauge theories. More precisely, strongly 
coupled QED in the noncompact lattice regularization 
has received considerable attention in recent time ( [^|J^] 
and references therein). Many interesting and exciting 



new results in the field of four-dimensional gauge theo- 
ries, as for instance the possibility to have a non Gaus- 
sian fixed point 0, appear in this model. The origin 
of these new phenomena is unfortunately not clear at 
all. It has been argued in the past that the monopole 
percolation transition could drive the chiral transition of 
strongly coupled QED Q. This issue is however rather 
unclear and it is therefore worthwhile to explore other 
possibilities. There are furthermore other exciting re- 
sults in non compact QED as the existence of a non 
chiral phase transition in the three-dimensional model 
H or the fact that the transition seems to prolongate 
to non vanishing fermion masses j^] which allow one to 
speculate that frustration could be the origin of these 
phenomena. We should say from the beginning that at 
the moment this is only speculation. Our aim is to in- 
vestigate this issue and, as a first step in this direction, 
we will study in this paper a model (jj in which frustra- 
tion is introduced by hand, although inspiration is taken 
from the hopping-parameter expansion of QED. Section 
is devoted to review how frustration appears in lattice 
QED. We analyze this phenomenon from the hopping pa- 
rameter expansion as well as from the loop expansion of 
the determinant of the Dirac-Kogut-Susskind operator. 
In Sec. ID we define our toy model, discuss its classical 
structure of ground states, introduce an order parameter 
for the phase transition and show how our model in two 
dimensions is equivalent to the XY-spin model. Section 
[V contains the numerical results of a simulation and the 
phase structure. The conclusions are reported in Sec. 



II. FRUSTRATION IN LATTICE QED 

We start with the action functional for QED on the 
lattice using Kogut-Susskind fermions, S = (3Sg + Sp, 
with Sg the pure gauge action (either compact or non- 
compact) and Sp the fermionic action, 
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V^(x)x(x){U^{x)x(x + fx) - U*(x - fi)x{x - n)} 



(1) 



where x labels the points and /i the forward unit vectors 
in the lattice. The E/ M (x) are compact U(l) link variables, 
arc fermion fields (Grassmann variables) and T}p,(x) 
are the Kogut-Susskind phases, defined as usual, 



Vnfr) = (-1) 



XX+X2 + ...+X,,- 



(1 < fi < d) 



(2) 



The fermionic part of the action is a bilinear of the 
fermion fields, Sp — x^x, with A = ml + A, and A 
an antihcrmitian matrix 

A *.y = \ ^2vn(x)[$x+n,yU f t(x) - 6 x - MtV U*(x - fj.)] (3) 

1 < M < d 

Due to its bilinear nature, we can integrate analytically 
the fermionic part of the action in the partition function, 
thus obtaining 



[dx][d X ][dU] 



-/3S G -xA x 



= J [dU] det A e~ pSG = J [dU] e~ s =« 

with the effective action defined as 

S eff = f3S G - In det A = /3S G - tr In A 



(4) 



(5) 



Now we are going to consider the loop expansion of the 
effective action, as well as the expansion of the determi- 
nant of the Dirac operator. 



A. The hopping parameter expansion 

Let us examine first the expansion of the effective ac- 
tion B. We have 

trln A = trln(m/ + A) = trln(m(/ + m _1 A)) 

= ^lnm + y(-l) fc+1 ^tr(A fc ) (6) 
^— ' k 

k>l 

The elements of A have a natural interpretation as di- 
rected links on the lattice: due to the local character 
of the action, only elements corresponding to pairs of 
nearest-neighborhood sites are non zero, so we can as- 
sign to each non-zero element A XjV the unique directed 
link which goes from x to y. We can thus rewrite the 
previous traces as 



tr(A fe ) = 

X,fJ, 1 ,fJ, 2 ---fik 

A^+pj^ h/ifc-ljS+jUi-l h/J/e-i+Mfc^PlH hAUe>° ( V 



In this paragraph only, and for the sake of clarity, the ^ 
represent backward as well as forward unit vectors in the 
lattice. 

We can assign to each non- vanishing term of the trace 
of A fe a unique closed and oriented loop on the lattice 
constructed with k links. Since loops must be closed, it 
is obvious that k has to be even in order to give a non- 
zero contribution to Consequently we rewrite the 
trace as 
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tr In A = V In » 



k>l 



(8) 



ieg-2 



where Gj represents all the loops with j links in the class 
considered, and A/ is a notation for the oriented product 
of the elements of A around loop I. The net result is 
real, as for the contribution of each loop we have the 
conjugate contribution of the loop taken in the opposite 
direction. To examine more closely this expression, and 
specifically the sign of the different terms, we use the 
definition of A. For each backward link we get a minus 
sign. But for a closed loop, exactly half of the links are 
backward links, so for a loop with 2k links we get a factor 
of (— 1) . On the other hand we have to take into account 
the Kogut-Susskind phases. One simple way to do this 
is to absorb them into the links, with a trivial change 
of variables, W^{x) = rj^ix^U^x). The corresponding 
change induced in the plaquettes consists in multiplying 
each of them by — 1. The final result for the effective 
fermion action expressed in the new variables W^(x) is 



SF ff = -V\*m + J2(-l? 



-2k 



k>\ 



2k(2) 2k ^ 



TV, 



(9) 



where Wi means the directed product of the W variables 
along loop I. We can see from this expression that the 
effective action is frustrated in the sense that there is no 
single configuration that minimizes term by term the ex- 
pansion. Different terms are minimized by different con- 
figurations. It should be said that even if we have used 
Kogut-Susskind fermions in our derivation, the conclu- 
sion remains valid for Wilson fermions. 

It is interesting to note that the analogous derivation 
for a scalar theory coupled to an Abelian gauge field does 
not show any kind of frustration. The essential differ- 
ences with respect to the fermionic case are on one hand 
that the corresponding determinant appears in the de- 
nominator, and on the other hand that in the scalar case 
the Kogut-Susskind phases are absent. The final result 
of the development of the effective action for the scalar 
case is the following: 



S ef f=(3S G + Vhim-^2 



-2k 



E u > ( 10 ) 



- 2k(2) 2k 
k>i y ' ieg 2k 

In this case each term of the expansion is minimized by 
the trivial configuration, so we have no frustration. 
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B. The expansion of the determinant 

The expansion discussed in the previous section is the 
most convenient approach for massive fermions. However 
one can always proceed directly with the expansion of the 
determinant of A, as long as the volume is kept finite 
So let us consider 

det A = det(mJ + A) = ^ £ lu i 2 ,---,i v A Ml • • • Ay iiv 

(11) 

where ^i 1 ,i 2 ,--,iv ^ s the totally antisymmetric symbol in 
V indices. Due to gauge invariance and to the proper- 
ties of the determinant, each non- vanishing contribution 
to can be associated to a graph on the L d lattice 
made up of some number of isolated sites (each one con- 
tributing with a factor m) and some number of closed, 
directed, simple (i.e., no link can appears twice) and non- 
intersecting loops. In other words, after eliminating the 
isolated sites, all the rest of the lattice should be covered 
with closed loops in such a way that in each site there are 
exactly one incoming and one outcoming link. In partic- 
ular this implies that the number of isolated sites (and 
so the powers of m) should be even. 

To minimize the contribution of the fermionic part to 
the effective action (||) one should look for configurations 
that maximize the fermionic determinant. In general, it 
is impossible to maximize simultaneously with a single 
configuration all the terms in the expansion of the deter- 
minant (this can be seen easily examining simple graphs 
containing two isolated sites, for example). 

III. THE TOY MODEL 

Motivated by the expansions analyzed in the previous 
section, we have introduced an Abelian gauge toy model 
with explicit frustration [Q. We have studied the effect 
of frustration in this model as a first step towards the 
analysis of the effect of frustration in more realistic mod- 
els, like QED. The model we have studied is a U(l) pure 
gauge model defined in a hypercubic lattice by the fol- 
lowing action: 

S = S 4 + S 6 = -0 RelJ pi + ReU P« ( 12 ) 

pi P6 

with (3, (3q > 0. The first part (S4) is the standard com- 
pact U{\) pure gauge Wilson action. The second part 
(Se) is defined as the sum of contributions over all closed 
loops made up with six non-repeated links. These loops 
can be generally classified in three different classes: pla- 
nar loops, loops that involve two planes and loops that 
involve three planes (obviously in the special case d = 2 
only planar loops appear). 



A. Classical ground states and order parameter 

It is not difficult to realize the existence of frustration 
in this model at the classical level in d > 2, even if we 
consider only the Sq piece of action (|l2]). In fact, it is not 
possible to minimize simultaneously the contributions to 
the action of all the loops. One could try to work out 
the generic minima for some of the classes of loops but 
instead, we will describe here, motivated by our numeri- 
cal results, a specific class of configurations that appear 
to be the classical ground states accessible to the system. 
This will permit us to define a natural order parameter 
for this system. 

In the configurations considered the plaquettes take 
the values ±1. It is easy to see that we can find configu- 
rations of this type which are minima for the planar and 
three-plane loops simultaneously. For the planar loops, 
the minimum condition requires that the sign of the pla- 
quette should be alternate on every plane (for this rea- 
son we will call these configurations antiferromagnetic, in 
analogy with magnetic systems). Among these configu- 
rations we can find a subset that minimize also the three- 
plane loops. There are a total of 2^ d_1 ^ configurations of 
this type, related by simple symmetry transformations. 
Now we should consider the two-plane loops. It is easy 
to check that these loops are not simultaneously mini- 
mized in the configurations that we have described. In 
fact, the net contribution to the action of these loops is 
zero in any of the minima considered. In the case d = 2 
only the planar loops appear in the action; therefore the 
antiferromagnetic configurations are global minima and 
the system is not frustrated. 

On the other hand the minimum for the S4 part of 
the action is reached for the homogeneous configuration 
U p i = 1, which becomes another source of frustration 
when considering the complete system. 

This analysis suggests the introduction of an order pa- 
rameter for each set of parallel planes, the staggered pla- 
quette, defined as follows: 

= ^E^) Rel W x )' e W = (-i) xi+x2+ - + - 

X 

(13) 

This order parameter is different from zero in the anti- 
ferromagnetic vacuum and vanishes in the ferromagnetic 
(homogeneous) one. For the configurations described be- 
fore it takes the values ±1. 



B. The case d = 2 and the XY model 

We will show in this subsection how, in d = 2, the 
thermodynamic limit of our model is the XY model. 

To this end let us consider first the Sq part of action 
(II2I). In two dimensions it contains only planar loops. It 
is a well known fact that in d = 2 it is possible to replace 
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the links by the elementary plaquettes as independent 
variables of the theory. This is due to the fact that there 
are no cubes in the lattice, and hence only one Bianchi 
identity: the product of all plaquettes must be equal to 
1. Therefore the number of independent plaquettes is 
V — 1. If we consider all the plaquettes as independent 
variables, the difference in the action is a term of order 
0(1), which is irrelevant in the large V limit. The pla- 
quette variables can now be assigned to the sites of the 
dual lattice, with an action which involves only nearest 
neighborhoods. Now, to obtain the usual form for the 
XY model, we still need another change of variables. We 
consider the usual bipartition of the planar lattice, and 
we change all the plaquettes of one of the sublattices to 
minus their complex conjugate. If we do this, the final 
form for the action is 



S 6 = p 6 J2 Re(Ui(-U*)) = -/3 6 ' 

<i,j> <i,j> 



(14) 



where X)<ij> means a sum over all pairs of nearest 
neighborhood sites of the (dual) lattice and 6k is the an- 
gle corresponding to the element of U(l) in each site. 
But this is just the action of the XY model on a square 
lattice. As is well known, this model has only one ground 
state and does not present frustration. 



The contribution of the S4 term to (12), expressed in 
the new variables, can be written as 



S 4 = -/?^Re(Ui) = ~/3^e iC os( 



(15) 



where equals ±1 according to which sublattice the site 
i belongs to. We see that the full action for the two- 
dimensional model, 



cos 



»)-/*E 



£i COS 



i) (16) 



<*,j> 



is the action of the X-Y model in the presence of an 
external staggered magnetic field. The X-Y model with 
an external random field has been extensively studied ( 
[ fi(|[i"l| and references therein) as a simple example of 
frustrated system, but we are not aware of studies of 
the staggered case, so we do not know to what extent 
frustration is relevant for this system. 



IV. NUMERICAL SIMULATIONS 

We performed Montecarlo simulations of this model in 
d = 3 and d = 4 using a standard Metropolis algorithm. 
The observables which were measured in the simulations 
were the staggered plaquette and the usual normalized 
plaquette in each plane, and the normalized contribution 
of the 6-loop part of the action, 



P ^ = ^ E <z)ReIV(x) , e(x) = ( _ ir +x 2 +-+x d 

1 < M < v < d 

^v = ^E Rel v ( 18 ) 

X 

p ^ = ^vJ2 RelJ ^ ( 19 ) 



where J\f means the number of 6-loops per volume. We 
also measured for completeness the imaginary part of 
these quantities, that gave results consistent with zero 
in all our simulations. 

In order to obtain a qualitative understanding of the 
phase diagram of this model, our strategy was to explore 
the parameter space over several lines of constant (3. To 
locate the possible transition points we did a number of 
annealing cycles over @§, while maintaining (3 fixed to a 
given value. The most extensive simulations were done 
in d = 4, and we discuss this case in more detail below. 




FIG. 1. Plaquette and 6-loop hysteresis cycles at (3 = 
against (3e 

As we have shown before, in d = 2 our model is equiv- 
alent to the X-Y model. It is well known that this model 
presents a continuous transition at a value of the coupling 
constant (3$ ~ 0.9, at least for (3 = 0. 

In Fig. |l| and | we show the results for the annealing 
cycles in d = 3 at (3 = 0. The figures suggest the exis- 
tence of a transition into a phase with a non-vanishing 
value of the staggered plaquette at a smaller value of the 
coupling constant than in d = 2, (3§ ~ 0.38. However, a 
more careful study and larger statistics are necessary to 
establish unambiguously the order of the transition. 



A. The system in d = 4 

We simulated the four-dimensional system along sev- 
eral lines of constant (3>0. First, we did several runs at 
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Fig. 2 

FIG. 2. Staggered Plaquette hysteresis cycle at f3 
against /3e 



Fig. 4 

= FIG. 4. Staggered Plaquette hysteresis cycle at /3 = 

against /fe 




Fig. 3 

FIG. 3. Plaquette and 6-loop hysteresis cycles at j3 = 
against /?6 



different (3, while maintaining /3 6 fixed to zero, in order to 
obtain properly thermalized configurations that we used 
later to start the annealing cycles in (3q. For each value of 
(3, we started the annealing with the thermalized configu- 
ration obtained in this way. After we finished the simula- 
tion for a given value of (3q, we changed it by some small 
amount and repeated the process starting with the last 
configuration generated at the previous value of (3§. Pro- 
ceeding in this way we did a complete forward-backward 
cycle. These simulations were done in 4 4 and 6 4 lattices. 

We show in Fig. ||| and ^ the resulting annealing cy- 
cles for some of the magnitudes we have measured at 
(3 = O.The staggered plaquette is shown here only for 
one of the planes (see the discussion below). We see a 
very clear hysteresis effect signaling a strong first-order 
phase transition. In order to check that this is not the 



effect of a bad thermalization around a continuous transi- 
tion we compared cycles obtained by increasing the sim- 
ulation time over two orders of magnitude, from 1000 to 
100000 Montecarlo sweeps, and no change in the width of 
the hysteresis region was seen. Also no significant finite- 
size effects were found in our simulations. We also im- 
plemented an over-relaxation procedure which was com- 
bined with the Metropolis algorithm; again, no change 
in the measured magnitudes was observed. Concerning 
the nature of the large (3q phase, we can see in figure 
4 that the staggered plaquette shows a sudden increase 
from zero, when entering into this phase. If we compare 
its value for the different planes we see that it only dif- 
fers in sign, with the signs being consistent with the ones 
obtained from the analysis of the c lassical ground state 
configurations done in section ( [II A ) . The simulations at 
larger values of (3q, produced configurations which were 
essentially the classical ground states described before. 

Now let us discuss our results for the lines at constant 
(3 > 0. For small values of f3, the results of the sim- 
ulations are qualitatively the same as for (3 = 0, apart 
from the fact that the transition point occurs at higher 
values of (3q for higher values of (3. But we know that 
this system, for /?g = 0, has a first-order transition at 
(3 ~ 1.0, the usual confined-deconfined transition for the 
standard compact U{\) pure gauge theory |]l^-|l4|. So it 
is very likely that a line of transition points emerges from 
the (3 axis and, indeed, this is what we found. We did 
simulations at several valued of (3 above the transition 
point. In Fig. |^ and ^ we show the results obtained by 
cycling over (3q while keeping (3 = 1.5. We can clearly see 
here two transitions: the first one should correspond to 
the continuation of the usual confined-deconfined transi- 
tion (the staggered plaquette remains zero for this transi- 
tion), whereas the second one corresponds to the transi- 
tion to the antiferromagnetic phase, as clearly shown by 
the jump of the staggered plaquette. The results for all 
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Fig. 5 

FIG. 5. Plaquette and 6-loop hysteresis cycles at f3 = 1.5 
against /?6 
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Fig. 6 

FIG. 6. Staggered Plaquette hysteresis cycle at (5 — 1.5 
against /3e 



the (3 values are consistent with this picture, and permit 
us to obtain a qualitative phase diagram, that we show 
in Fig. 0. We found three phases separated by two first- 
order lines. The antiferromagnetic phase is characterized 
by a non-zero value of the staggered plaquette, which is 
zero in the other two phases. In this phase we have sev- 
eral states, to be precise eight states, related by a spon- 
taneously broken translational symmetry. The other two 
phases are the continuation of the confined and uncon- 
fined phases of the standard compact pure gauge model. 



V. COMMENTS 

Our starting point was the observation that compet- 
ing interactions are present in the effective gauge theory 
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obtained, after integration on Grassmann fields, from a 
gauge-fermion theory on the lattice, such as QED. In- 
spired from analogies with condensed matter models like 
spin-glasses in which frustration and disorder conspire 
to give non standard phases to the system, we asked if 
the rich, and in some sense unexpected fenomenology ob- 
served in 3d and Ad lattice QED can be ascribed to the 
presence of non standard vacua. 

The first step in our analysis has been the study of 
a simpler model, in which frustration is introduced by 
hand, but as suggested by the hopping parameter expan- 
sion. In the d = 3 and d = A cases we have seen clear 
signals of the existence of an antiferromagnetic phase in 
which translational invariance is broken. This is not to- 
tally surprising since the solutions of classical equations 
of motion, at least in a limiting case of the parameter 
space, share this property. A crucial point in the analy- 
sis has been the availability of a clear and simple order 
parameter, inspired from the properties of these classical 
configurations. The d = 2 case was found analytically 
equivalent to the X — Y model with staggered external 
magnetic field in the thermodynamical limit. 

The simple model we have considered has only a finite 
number of competing interactions, so on general grounds 
we cannot expect to find spin glass-like phases. On the 
contrary, in the full model, the number of frustrated 
terms in the effective action is diverging with the vol- 
ume of the system. We do not know if this is enough to 
simulate in some way the inclusion of disorder, that in a 
spin glass system is the other ingredient for the appear- 
ance of non standard phases: this is an issue for future 
work. The first problem we have to face is the definition 
of a suitable order parameter. 

The list of features in three and four dimensional QED 
that lack an explanation is long, and we suspect that, at 
least some of them, could be related to frustration at the 
level of the effective action. 
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